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Abstract 

This paper is devoted to consideration of the theory of collisionless 
statistical systems with interparticle scalar interaction. The mathemati¬ 
cal model of such systems is constructed and the exact solution of Vlasov 
equation for isotropic homogenous model of the Universe is found. Asymp¬ 
totic solutions of self-consistent Vlasov - Einstein model for conformally 
invariant scalar interactions are found. 
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1 Introduction 

In series of Author’s works n, 0, m etc. the mathematical model of self- 
gravitating statistical system of scalar charged particles is formulated. In par¬ 
ticular, in [3] this model is extended on sector of negative effective masses of 
scalar charged particles. However, in these papers as well as in other works 
of the Author with his pupils it was usually considered a locally equilibrium 
system. As is well-known, at inflationary stage thermodynamic equilibrium is 
violated [5] hence it is reasonable to consider collisionless situation realized at 
this stage which is described by corresponding Vlasov equations. This paper is 
devoted to this problem. 


2 Self-consistent Mathematical Model of Colli¬ 
sionless Plasma of Scalar Charged Particles 


2.1 Collisionless Kinetic and Transport Eqnations 


Distribution functions fa{x,P) of particles of “a” sort having scalar charge qa 
are defined by invariant collisionless kinetic equations m 


[HaJa]=0, 


( 1 ) 


where [Ha, fa] is a Poisson bracket which can be written using covariant Cartan 
operator of differentiation V in fornjl: 


[H,^] 


dH~ ^ ST - 

-V,T-Vjil, 

dP^ " SP. * ’ 


^i.e. Vlasov equations. 
^Details see in [4] 


( 2 ) 
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Hamilton function of a particle in scalar field and normalization ratio for gen¬ 
eralized momentum have the next form: 

Hix,P) = ^[m-^{x){P,P)-m^] =0, (3) 

iP,P)=ml = {qga^f. (4) 

Distribution functions fa{x,P) of particles of “a”-sort having scalar charge 

Qa are defined through invariant collisionless kinetic equations [1]: 

[i?a,/a]=0, (5) 


where Hamilton function and normalization ratio for generalized momentum 
have the next form: 

H{x,P) = ]^[m~^{x){P,P)-m^] =0, (6) 

{P,P)=ml = {q,a^f. (7) 

dH ~ ~ 

Let us notice identity laws, valid for Hamilton function ([6]) which will be 
useful in future: 

VjiL = -Vim*, (9) 

= + ( 10 ) 

//fc* Ui i 

where 4'(a;,P) is an arbitrary function. 

Further, let us define macroscopic moments using distribution functions: 

29-1-1 f 

na = 9a J fa{x,P)P"dPo. (11) 

Po 

- particle number density vectors which are identically conserved in consequence 
of kinetic equations ([5|) ; plasma energy-momentum tensor 

= E P)P^P'^dPo (12) 

“ Po 

and plasma scalar charge density 

= ^E ^(2^)3^ ^° / fa{x,P)dPo. (13) 

“ Po 

These macroscopic moments are related to each other by transport equations 
(energy-momentum conservation law). 


VfcT;'= - crV*$ = 0. 


(14) 


2 





2.2 Scalar Field Equation 

First, let us consider Lagrangian function for a classical massive real scalar field 
$. In such a case the Lagrangian of real scalar field can be chosen in the next 
form: 

Ls = ^ , ( 15 ) 

where rus is a mass of scalar field’s quanta and it is £2 = 1 for a classical scalar 
field, £2 = ~1 for a phantom scalar field; £1 = 1 for a field with repulsion of like- 
charged particles, £1 = —1 for a field with attraction of like-charged particles 
£1 = — 1. Then energy-momentum tensor of a scalar field is: 

rf = |1 ( 2 $’*$’'= - j -p £2g*'=m2$2) . (16) 

Let us write the equation of massive non-conformal scalar field with a source 

[U: 

□d)-I-mg<i) = —47r£i(T, (17) 

where 

□d> = ^ 

\/—g ox^ ox'^ 

Let us consider now the Lagrangian function for a classical massive real 
conformal scalar field $ (see e.g., [4,5]; for a massive scalar field the conformal 
invariance is understood as asymptotic property at (rn-s —>-0)): 

is = - £2 w2$2 + , (18) 

The Lagrangian function differs from the standard one (see e.g. i) in presence 
of factor I/Stt as well as introduced unit indicators £q,. In addition, Ricci tensor 
in the article is obtained by means of convolution of first and third indices of 
Riemann tensor Rji = Components of energy-momentum tensor of a 

scalar field relative to Lagrangian function (IT^ are equal to [HIT]: 


rpik _ 

" “ Stt 


+ £2771^3*'"$^ -I- 

3 3 


- ]:Rg"’")^‘^ - -f 

o z 00 


(19) 


assuming $ ^ Const, we get equation of massive scalar field with a source (see 





2.3 Einstein Equations 

Complete system of macroscopic equations comprises of kinetic equations (O, 
field equations (HZl) or (EOI) and Einstein equations: 

R"'" - \Rg"’" = 87r(r;'= + Tf), (21) 

where is a defined above energy-momentum tensor of a statistical system 
and is an energy-momentum tensor of a scalar field dm) or dm). Covari¬ 
ant divergence of Einstein equations turns them into identities as a result of 
transport equation da and field equation da or da- 

3 Cosmological Evolution Equations 

3.1 Kinetic Equations 

Let us assume scalar potential to be function only on time $(??) for a space-flat 
Friedmann Universe: 

ds^ = — dx^ — dy"^ — dz^) = a?'{vi)dsQ. ( 22 ) 

Let us also assume distribution function depending only on time and squared 
space momentum: 


fa = fa{V,P); +P|+P|, (23) 

so that kinetic energy of particles is equal to 

E = = — = V + 3- (24) 

a \ 


Calculating Poisson bracket da, let us reduce kinetic equations dSD to the next 
form: 



(25) 


Thus, homogenous isotropic solution of collisionless kinetic equation in Fried¬ 
mann metrics does not explicitly depend on time variable which was also true 
for a gas of neutral particles mm- 


/(r;,Po)=/(Po). 


(26) 


3.2 Sewing of the Collisionless Solution and Equilibrium 
Solution 

To find certain distribution in plasma where previously thermodynamic equilib¬ 
rium was previously maintained, let us apply next method [319]. Let us assume 
that till certain instant rjo local thermodynamic equilibrium was maintained in 
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plasma and it was instantly violated at the same instant. Thus, at rj < rjQ 
distribution functions were locally equilibrium i.e.: 


+ 


foiv^Pi) = 

a9{r]) 


exp^-- 
n -1 

± 


Kv) , 
e 

V < Vo 


(27) 


(/r are chemical potentials, 0 is a temperature), and at rj >= rjo they were 
described by collisionless kinetic equations. Therefore at instant v = Vo '■ 


fiPo) 



\/g^a^(bo)^^(bo) + PqI 

Sivo)a{vo) 


± 1 


Kvo) 

S{vo) 


V = Vo- 


(28) 


The single option to fulfill this condition is possible if Pq in equilibrium distri¬ 
bution is substituted by integral: 


lp2(77,Po)=Po' 



a — — drj. 
dr] 


Vo 


(29) 


Thus, the exact solution of collisionless kinetic equations (IS|), which becomes 
equilibirum one at rj = tjq, is: 


f{v,Po) = 


f Mo , 

Vmo+P'^\ 1 / 


00 ) \ 


1 -1 


(30) 


where the following denotations are introduced: 


Too = aoq^o; ao = a(ryo); Mo = Kvo); 
do = d{vo)] ^"0 = ^(bo); Po = aoP- 


(31) 

(32) 


Thus, mo is an effective particle mass at instant of equilibrium violation. Further 
for the sake of brevity we will call the solution (1301) quasi-equilbrium solution. 


3.3 Moments of Quasi-equilibrium Distribution Function 

It is easily seen that quasi-equilibrium distribution function’s moments mD 
and m which are interesting for us, coincide in the algebraic structure with 
corresponding moments of locally-equilibrium distribution: 


n* = u^n] u* = -51] 

(33) 

a 

= {Spi + Ppi)u^u'^ - Ppig^\ 

(34) 


where u® is a unit timelike vector of a synchronous observer. However, macro¬ 
scopic scalars n, P, £, a relative to quasi-equilibrium distribution ((3(111 do not 
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coincide with corresponding scalars relative to equilibrium distribution func¬ 
tion. Calculating, we fined: 


pal 


p^dp 


2TT^a^{ri) 


exp, -EM 


So 


± 1 


(35) 


so that: 


na^ = Const, 


(36) 


which ensures conservation law of particle number. Thus, the expression for 
particle number density coincides with standard equilibrium one with constant 
mass, temperature and chemical potential within the accuracy of multiplication 
by scale factor (ao/a)^. 

Expressions for remaining three scalars do not coincide with equilibrium 
ones: 


£pi = 



p^^Jm? -I- p'^dp 

Km) _|_ ^ ^ 

^0 ^0 j 


(37) 


Ppi = 



dp 


+P" expf-4^ + ± 1 ’ 


( \ 2 2 
ao\ pq^ 

-) 

dp 


y/m2+p2exp|^_^ ^ ± 1 ’ 


(38) 


(39) 


where: 

m = rriA ,— = —^rno = m[9). (40) 

oo ao'Po 

The behavior of these scalars is substantially defined by the behavior of invari¬ 
ant: 

$ = a{p)^{p). (41) 


Let us note that at 4) = Const integrals in expressions for all three scalars are 
constant values. 

= p,/9 which is a reduced chemical potential 
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3.4 Field Equations 

In case of conformally non-invariant field m field equation (I17|) in Friedmann 
metrics takes form: 


-I-£27715$ =-dTreicr. (42) 

dr] dr] 

In case of conformally invariant scalar field (I18p this equation needs to be sub¬ 
stituted by the following one: 


-I- £27715$ = -47r£itT. 
dr]-^ 


(43) 


Let us consider conformally invariant case of massless scalar field. Let us 
rewrite scalar charge density (I39|) in the following form: 


cr = 


cr($) 


(44) 


where 






p 


dp 


-C V^o+P^ 

1 00 Bo 


± 1 


(45) 


Thus, conformally invariant equation (1461) for a massless scalar field takes form: 


-—$ = —47r£l(T($). 
dr]^ 


(46) 


i.e. equation for scalar conformally invariant massless field with a source be¬ 
comes autonomous ordinary second-order integro-differential equation. We can 
find the solutions of this equation in different extreme cases. 

Let us first consider a ultrarelativistic case: 


(p) ^ max(mo,m). 


(47) 


In such a case equation (HSl) becomes a linear second-order differential equation 
with constant coefficients: 


where 


^ - 

dp' 


2 $ -k £10;^$ = 0, 


f pdp 

27 r 2 J exp{ -^;°+P ) ± 1 ■ 


(48) 


(49) 
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In a case of like-charged particles repulsion (ei = -1-1) this equation has its 
general solution: 


^ ^ cos(a;o77) , ^ sin(a;o?7) 

^ — <-"i--f —, 

a[V) a(r]) 

and in a case of attraction (ei = — 1 ) it is: 

^ ^ cosh(wo?7) , ^ sinh(a;o7?) 

f — Cl- - -h 02- l-T — 

a[r]) a{r]) 


(50) 


(51) 


Let us notice that exactly the same solution was found in m for ultrarelativistic 
equilibrium plasma. This coincidence is easily explained since quasi-equilibrium 
distribution (|30|) coincides with equilibrium one in the ultrarelativistic limit. 
Let us consider now a non-relativistic cas^: 


{p) <C m. (52) 

Then instead of homogenous equation (|4^ we get inhomogeneous equation: 


where 


CTO = 


drj 




= _47reicro, 


pdp 


27r^mo J „ ^ / —muo + y/m^+p'^ ^ 

0 exp(^ ^ 

The solutions of this equation are: 


a a a 


(53) 


(54) 


(55) 


3.5 Einstein Equation 

It can easily be seen that energy-momentum tensor of conformally invariant 
massless scalar held $( 77 ) has a form of energy-momentum tensor of the ideal 


hux (ei = -|-1) 

= {£s + PsVv^ - Psg^\ 

(56) 

where 


(57) 


- 24^a4 + 8^ei$a($)] . 

(58) 

so that: 

Ts = £s- 3P, = 

(59) 


^Let US note that here the smallness of average momentum is required relative to m but 
not relative to mo- 













The single non-trivial Einstein equation in Friedmann metrics takes the next 
form: 

^ = ^i£,+£^,). (60) 

Qi o 

Substituting the expressions for energy density from (1571) and dSTD into this 
equation, we reduce it to: 




Stt 








- p^dp 


0 exp 


^0 


00 


± 1 


(61) 

(62) 


Thus, cosmological model in case of massless conformally invariant scalar field 
with a collisionless source is reduced to two ordinary integro-differential equa¬ 
tions of the second order (H51) and first order (EH) relative to functions 0 ( 77 ) and 
#( 77 ). 

In particular, for the ultrarelativistic case we have £pi(^) = £q = Const. 
Choosing a convergent solution of scalar field equation (IbH]) 


$ = ())o sin(a;o7?) 


(63) 


and substituting it into Einstein equation (1611) . we reduce it to the next form 


/2 Stt ,99 2 r* 

a = —(PqLUqCOS 07077 -I-—C 0 . 

The solution of this ordinary differential equation has the following form: 


(64) 


fc" = 


a{p) = ^ 


(65) 


2 2 /q 2 c 

a = 00 ^ 0 ; P = '^ 0 , 

( 66 ) 

q,2 -(_ ^2 ’ 


E(a;, k) is an elliptic integral of the second kind: 

X 

E(a 7 o 77 , “ y \/1 — sin^ xdx. 


(67) 
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